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Abstract 

The novel point of this analysis is a direct use of reggeon exchange technique for the 
description of the reactions ir N — > two mesons + N at large energies of the initial pion. 
This approach allows us to describe simultaneously distributions over M (invariant mass 
of two mesons) and t (momentum transfer squared to nucleons). Making use of this 
technique, the following resonances (as well as corresponding bare states), produced in 
the ttN -> vr 7r°iV reaction are studied: / (980), /o(1300) (/o(1370) in PDG notation), 
/o(1200 - 1600), / (1500), /o(1750), / 2 (1270), / 2 (1525), / 2 (1565), / 2 (2020), / 4 (2025). 
Adding data for the reactions pp(atiest, fromliquidH 2 ) — ► 7r 7r°7r°, 7r°7r°?7, Tr°r]r] and 
pp(atrest, from gaseous H 2 ) — ► 7r°7r 7r°, 7r°7r°r/, 7r°ryry, we have performed simultaneous K- 
matrix fit of two- meson spectra in all these reactions. The results of combined fits to 
the above-listed isoscalar /j-states and to isovector ones, ao(980), ao(1475), a 2 (1320), are 
presented. 

PACS numbers: 11.25.Hf, 123. IK 



1 Introduction 

The study of the mass spectrum of hadrons and their properties is the key point for the 
understanding of colour particle interactions at large distances. But even the meson sector, 
though less complicated than the baryon one, is far from being completely understood. We 
mean that 

(i) there is no sufficient information about states above 2 GeV, 

(m) certain quark-ant iquark states below 2 GeV (e.g. 2 states) are still missing, 
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(Hi) there is no clear understanding of the glueball spectrum (although strong candidates in 
the ++ and 2 ++ sectors exist, we have no definite information about the h sector), 

(iv) some analyses reported the observation of other exotics (e.g. hybrid) states, 

(v) in the scalar sector not only the properties but also the existence of states like a, k, /o(1300) 
(/o(1370) in PDG notation) is under discussion. 

So, there is indeed a strong demand for new data which can help us to identify the meson 
states in a more definite way. However, the situation is only partly connected with the lack 
of data. In the lower mass region there is a lot of data taken from the proton-antiproton 
annihilation at rest (Crystal Barrel, Obelix), from the 77 interaction (L3), from the proton- 
proton central collisions (WA102), from J/ty decay (Mark III, BES), from D- and 5-meson 
decays (Focus, DO, BaBar, Belle, Cleo C) and from ttN — > two mesons + N reactions with 
high energy pion beams (GAMS, VES, E852). Most of these data are of high statistics, thus 
allowing us to determine resonance properties with a high accuracy (though, let us emphasize, 
in the reactions ttN — > two mesons + N polarized-target data are lacking). 

Nevertheless, in many cases there are significant contradictions between analyses performed 
by different groups. The ambiguities originate from two circumstances. 

First, in the discussed sectors the analyses of data taken from a single experiment cannot 
provide us with a unique solution. A unique solution can be obtained only from the combined 
analysis of a large set of data taken in different experiments. 

Second, there are some simplifications inherent in many analyses. The unitarity was ne- 
glected frequently even when the amplitudes were close to the unitarity limit. A striking 
example is that up to now there is no proper i^-matrix parametrization of the 1 and 2 ++ 
waves which are considered by many physicists as mostly understood ones. As to multipartical 
final states, only a few analyses have ever considered the contributions of triangle or box singu- 
larities to the measured cross sections. However, these contributions can simulate the resonant 
behavior of the studied distributions, especially in the threshold region (for more detail, see [I] 
and references therein). 

In the analysis of meson spectra in high energy reactions ttN — > two mesons + N, many 
results are related to the decomposition of the cross sections into natural and unnatural ampli- 
tudes that is based on certain models developed for the two-pion production at small momenta 
transferred, (e.g., see [21 El II])- However, as was discussed by the cited authors, a direct applica- 
tion of these methods at large momenta transferred to the analysis of data may lead to a wrong 
result. In addition, the ttN — * two mesons + N data were discussed mostly in terms of t-channel 
particle exchange, though without proper analysis of the t-channel exchange amplitudes. 

A decade ago our group performed a combined analysis of data on proton-antiproton anni- 
hilation at rest into three pseudoscalar mesons, together with the data on two-meson S"-waves 
extracted form the ttN — > ttttN, nnN, KK and nn'N reactions [5j El [7J. The analysis has 
been carried out in the framework of the i^-matrix approach which preserves unitarity and 
analyticity of the amplitude in the two-meson physical region. Although the two-meson data 
extracted from the reaction irN — > two mesons + N at small momentum transfer appeared to 
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be highly compatible with those found in proton-antiproton annihilation, we have faced a set 
of problems, describing the irN — > two mesons + N data at large momentum transfer. As we 
have seen now, the problems were owing to the use of partial wave decomposition which was 
performed by the E852 Collaboration and showed a huge signal at 1300 MeV in the S'-wave. 

The strategy of our present approach is as follows. The analysis of a large set of experimental 
data on proton-antiproton annihilation at rest is carried out together with the analysis of the 
irN — > two mesons + N data based on the t-channel reggeized exchanges. For the irN — > 
two mesons + N reactions, the data at small and large momentum transfers are included. Here, 
as the first step, we perform the analysis in the framework of the .ff -matrix parametrization 
for all fitting channels (.ff-matrix approach insures the unitarity and analyticity in the physical 
region). At the next stage, we plan to use the N/D method for two- meson amplitudes satisfying 
these requirements in the whole complex plane. 

In this paper, we present the method for the analysis of the irN interactions based on the 
t-channel reggeized exchanges supplemented by a study of the proton-antiproton annihilation 
at rest. The method is applied to a combined analysis of the nN — > tt tt N data taken by E852 
at small and large momentum transfers and Crystal Barrel data on the proton-antiproton 
annihilation at rest into three neutral pseudoscalar mesons. The even waves, which contributed 
to this set of data, are parametrized within the i^-matrix approach. To check a strong S'-wave 
signal around 1300 MeV, which has been reported by E852 Collaboration from the analysis of 
data at large momentum transfers, is a subject of a particular interest in the present analysis. 

We present the results of the new i^-matrix analysis of two-meson spectra in the scalar, 
J p = + , and tensor, J p = 2 + , sectors: these sectors need a particular attention because just 
here we meet with the low-lying glueballs, / (1200 — 1600) and / 2 (2000). The situation with 
the tensor glueball is rather transparent allowing us to make a definite conclusion about the 
gluonium structure of ^(2000), while the status of the broad state /o(1200 — 1600) requires a 
special discussion: this state is nearly flavour-blind but the corresponding pole of the amplitude 
dives deeply into the complex-M plane. It is definitely seen only in the analysis of a large number 
of different reactions in broad intervals of mass spectra (for example, see [T] and references 
therein) . 

So, here we consider the following reactions: 

(i) irp — > TT°TT°n at high energies of initial pion and small and large momenta transferred to 
nucleon, and 

(ii) pp(atrestH2) — > 7r 7r°7r°, TT°7T°n, ir°nn in liquid and gaseous H2 — the data on these reac- 
tions give us the most reliable information about scalar and tensor sectors. 

As was stressed above, the novel point of the performed .fT-matrix analysis is the use of reggeon 
exchange technique for the description of irp — > ixixn at high energies that allows us to analyze 
the two-meson invariant mass spectra and nucleon momentum transfer distributions simulta- 
neously. 

The paper is organized as follows. 

In Section 2, we consider meson-nucleon collisions at high energies and present formulas 
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for peripheral two-meson production amplitudes in terms of reggeon exchanges. Amplitudes 
for the description of low-energy three-meson production in the i^-matrix approach are given 
in Section 3. The fitting procedure is described in Section 4. In Conclusion we summarize the 
results. Technical aspects of the fitting procedure are discussed in [8]. 



2 Meson— Nucleon Collisions at High Energies: Periph- 
eral Two-Meson Production in Terms of Reggeon Ex- 
changes 

The two-meson production reactions up — > -Kirn, KKn, nnn, nn'n at high energies and small 
momentum transfers to the nucleon are used for obtaining the S"-wave amplitudes tttt — > tttt, 
KK, rjn, nn' at |t| < 0.2 (GeV/c) 2 because, as commonly believed, the tt exchange dominates 
this wave at such momentum transferred. At larger momentum transfers, \t\ > 0.2 (GeV/c) 2 , 
we observe definitely a change of the regime in the S-wave production — a significant con- 
tribution of other reggeons is possible (ai-exchange, daughter- tt and daughter-ai exchanges). 
Nevertheless, the study of the two-meson production processes at \t\ ~ 0.5 — 1.5 (GeV/c) 2 looks 
promising, for at such momentum transfers the contribution of the broad resonance (the scalar 
glueball /o(1200 — 1600)) vanishes. Therefore, the production of other resonances (such as the 
/o(980) and /o(1300)) appears practically without background - this is important for finding 
out their characteristics as well as a mechanism of their production. 

What we know about the reactions irp — > -Kirn, KKn, nnn, rjrj'n allows us to suggest that a 
consistent analysis of the peripheral two-meson production in terms of reggeon exchanges may 
be a good tool for studying meson resonances. Note that investigation of two-meson scattering 
amplitudes by means of the reggeon exchange expansion of the peripheral two-meson production 
amplitudes was proposed long ago [H] but was not used because of the lack of data until now. 

The i^-matrix amplitude of the peripheral production of two mesons with total angular 
momentum J reads: 

^N{k3)G R 4j N {p 2 )^jR{s 7TN ,t)K 7Tm {s) l-p{s)K{s) ' Q {J) {k u k 2 ) , (1) 

This formula is illustrated by Fig. [1] for the production of tttt, KK, nn, nn' systems. Here the 
factor (ipNi^GRipN^)) stands for the reggeon-nucleon vertex, and Gr is the spin operator; 
R(s n N, t) is the reggeon propagator depending on the total energy squared of colliding particles, 
SnN = G°i +P2) 2 , and the momentum transfer squared t = (p2 — k 3 ) 2 , while the factor K nR M [1 — 
ip(s)K(s)]~ 1 is related to the block of two-meson production; s = M 2 = (ki + k 2 ) 2 , and 
p(s) is the phase space matrix . In the reactions up — > irnn, KKn, nnn, nn'n, the factor 
K n R(t)(s)[l — zp(s)i^(s)] _1 describes transitions irR(t) — > hit, KK, nn, nrj': in this way the 
block K n R( t ) is associated with the prompt meson production, and [1 — ip(s)K(s)]~ 1 is the 
K-matrix factor for meson rescattering (of the type of 7T7t — > tttt, tttt — > KK, KK —>■ nn, and 
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Figure 1: Example of a reaction with the production of two mesons (here tttt, KK, nn, nn' in 
ir~p collision) due to reggeon [K] exchange. 

so on). The prompt-production block for transition ttR — > b (where b = tttt, KK, nn, nn', An, 
...) is parameterized with singular (pole) and smooth terms [5| [71 [TO]: 

SUo).„ = E^#*- + WM)- (2) 

The pole singular term, G^(t)g^ / (/i^ — s), determines the bare state: here G^(t) is the bare 
state production vertex while the parameters and /i n are the coupling and the mass of the 
bare state - they are the same as in the partial wave transition amplitudes nit — > nn, KK, 
rjrj, rjrj 1 , An, .... The smooth term f n B.,b stands for the background production of mesons. The 
G%x(t), f n R y b, gi j are free parameters of the fitting procedure, while the characteristics of 
resonances are determined by poles of the .ff -matrix amplitude (remind that the position of 
poles is given by zeros of the amplitude denominator, det\l — p(s)K(s)\ = 0). 

Below we explain in detail the method of analysis of meson spectra using as an example the 
reactions ttN — > titiN , KKN , rjrjN, rjrj'N, nnnnN . 

2.1 Reggeon exchange technique and the K- matrix analysis of me- 
son spectra in the waves J PC = ++ , 1 , 2 ++ , 3 , 4 ++ in high 
energy reactions irN — > two mesons + N 

Here we present the technique of the analysis of high-energy reaction ir~p — > mesons + n, with 
the production of mesons in the J PC = ++ , 1 , 2 ++ , 3 , 4 ++ states at small and moderate 
momenta transferred to the nucleon. 

The following points are to be emphasized: 

(1) The technique can be used for performing the i^-matrix analysis not only for ++ and 2 ++ 
wave, as in [5l[7J[T0], but simultaneously in 1 , 3 , 4 ++ waves as well. 

(2) We use the reggeon exchange technique for the description of the t-dependence in all an- 
alyzed amplitudes. This allows us to perform a partial wave decomposition of the produced 
meson states directly on the basis of the measured cross sections without using the published 
moment expansions (which were done under some simplifying assumptions - it is discussed 
below in more detail). 

(3) The mass interval of the analyzed spectra is extended up to 2500 MeV thus overlapping 
with the mass region studied in reactions pp(in flight)-^ mesons [TT] . 



5 



We discuss in detail the reactions at incident pion momenta 20-50 GeV/c, such as measured 
in [121 H31 HH [151 HH H7] : (i) 7r~p -»• Tr+Tr^+n, (ii) n~p -> 7r°7r° + n, (Hi) n~p -> K s K s + n, (iv) 
ir~p — > rjr] + n. At these energies, the mesons in the states J PC = ++ , 1 , 2 ++ , 3 , 4 ++ are 
produced via t-channel exchange by reggeized mesons belonging to the leading and daughter 
7r, a,\ and p trajectories. 

But, first, let us present notations used below. 



2.1.1 Cross sections for the reactions nN — > mrN, KKN, r)T)N 

We consider the process of the Fig. [1] type, that is, irN interaction at large momenta of the 
incoming pion with the production of a two-meson system with a large momentum in the beam 
direction. This is a peripheral production of two mesons. 

The cross section is defined as follows: 

(2tt) 4: \A\ 2 

da = d(j)(px + p 2 ; fa, fa, fa), 

°V SttN\P2 I cm (np) 

d<P(pi + p 2 , fa, fa, fa) = (27r) 3 d$(P; fa, fa) d$( Pl + p 2 ; P, fa) ds , (3) 

where \p2\cm(-Kp) is the pion momentum in the cm. frame of the incoming hadrons. Taking into 
account that invariant variables s and t are inherent in the meson peripheral amplitude, we 
rewrite the phase space in a more convenient form: 

— * 

Momentum |/ci| cm (i2) is calculated in the cm. frame of the outgoing mesons: in this system one 
has P = (M, 0,0,0, ) = (a/s, 0,0,0) and g^„fa u = —g^fau = (0, k sin 6 sin 92, k cos 6 sin y?, 
k cosGfc) while dfl = d(cosQ)dip. We have: 

(2vr) 4 |A| 2 (2vr) 3 1 dt dM 2 d$(P, fa, fa) _ \A\ 2 p(M 2 ) MdM dt dQ 

8\P2\cm(np)y/s^N (2tt) 5 8\p 2 \ cm (np) y/s^ 32 ( 27r ) 3 |P2|L(7rp) 8 *N 

The cross section can be expressed in terms of the spherical functions: 

^ = N(M, t) Z (prWfrr) + 2± (YDReYne,,)). (6) 

The coefficients N(M,t), (Y®), (Y™) are subjects of study in the determination of meson 
resonances. 

Before describing the results of analysis based on the reggeon exchange technique, let us 
comment methods used in other approaches. 
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2.1.2 The CERN-Munich approach 



The CERN-Munich model \Tb\ was developed for the analysis of the data on the n~p — ► 
7T + 7T~n reaction. It is based partly on the absorption model but mainly on phenomenological 
observations. The amplitude squared is written as 



\A\ 



J=0 J=l J=l 



(7) 



and additional assumptions are made: 

1) The helicity-1 amplitudes are equal for natural and unnatural exchanges A^ = A^ 



2) The ratio of the A^ and the A® amplitudes is a polynomial over the mass of the two-pion 

f-1 n ( 3 X ~ 1 

system which does not depend on J up to the total normalization, Aj = AjlCj J2 b n M r 

\ n=0 

Then, in [15], the amplitude squared was rewritten using density matrices = A^A^, 
PS? = A° n Atf*, p?F = 2A^Atf as follows: 



d 



i,k,l 
n,m,J 



J=0 \n,rn / 

+ E ReY}(&, if) (E <° m >™ + d°n£jPn 

J=0 \n,m 

J dtt ReY*(Q, <p) ReY*(e, <p) ReYj(e, tp) 



J dQ ReYj(Q, tp) ReY}(0, <p) 



(8) 



Using this amplitude for the cross section, the fitting to the moments < Y™ > has been carried 
out. 

The CERN-Munich approach cannot be applied to large t, it does not work for many other 
final states either. 



2.1.3 GAMS, VES, and BNL approaches 



In GAMS p^2l [13] . VES [16], and BNL [T7] approaches, the ttN data are described by a sum of 
amplitudes squared with an angular dependence defined by spherical functions: 



\A 2 



E^i , (e,^) + E4~ ) v / 2i?eF 7 1 (e,^) + ^A™y/2ImY}{e,<p) 



(+). 



(9) 



The A j functions are denoted as So, Pq, D , F . . ., the Aj ^ functions are defined as P_, £)_, F_, . . 
and the A^ functions as P + , D + , F + , . . .. The equality of the helicity-1 amplitudes with natural 
and unnatural exchanges is not assumed in these approaches. 

However, the discussed approaches are not free from other assumptions like the coherence 
of some amplitudes or the dominance of the one-pion exchange. In reality the interference of 
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the amplitudes being determined by t-channel exchanges of different particles leads to a more 
complicated picture than that given by (jUJ), this latter may lead (especially at large t) to a 
misidentification of quantum numbers for the produced resonances. 

For example, in [17] the S-wave appears in an unnatural set of amplitudes only. Natural 
exchanges have moments with m=l,2,3.... However, the al-exchange is a natural one, therefore 
it contributes into the S-wave and does not interfere with unnatural exchanges - in this point 
the moment expansion [TTJ does not coincide with formula with reggeon exchanges. 



2.2 The ^-channel exchanges of pion trajectories in the reaction 

7T~p — > 7T7T Tl 

Consider now in more detail the production amplitude for the irir system with 1 = and 
J PC = ++ , 2 ++ and show the way of its generalization for higher J. 



2.2.1 Amplitude with leading and daughter pion trajectory exchanges 

The amplitude with t-channel pion trajectory exchanges can be written as follows: 

A^lTn toncs) = E aUr{^) -> vnrW^g 2 ) 0+(^1)</?p) 9%\t) (10) 

The summation is carried out over the leading and daughter trajectories. Here A(7rR(irj) — > tttt) 
is the transition amplitude for meson block in Fig. (TJ g^i is the reggeon-iViV coupling and 
RkAsttn, q 2 ) is the reggeon propagator: 



R^ 3 {s nN ,q )=exp(-t-a^(q )) , jz {j) , (l {j) ON — ■ (Hj 



■ 7T ia\, o,\ {SttN/ SttNO) 

sin (f o4% 2 )) T (Iq^V) + 1 

The 7r-reggeon has a positive signature, £ n = +1. Following [TJ [181 EH EQ], we use for pion 
trajectories: 

a^ eadin9 \q 2 ) ~ -0.015 + 0.72g 2 , a (d™ 9 hter-i) ^ _ _ 1AQ + g.72g 2 , (12) 

where the slope parameters are given in (GeV/c) -2 units. The normalization parameter s^ato 
is of the order of 2-20 GeV 2 . To eliminate the poles at q 2 < we introduce Gamma- functions 
in the reggeon propagators (recall that 1/T(x) = at x = 0, —1, —2, . . .). 

For the nucleon-reggeon vertex we use in the infinite momentum frame the two- 
component spinors <p p and ip n (see, for example, [T| [T8l I2T] ): 

gSihh&ip*)) — {<pt(<?<T±)<Pp) C\t) . (13) 
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As to the meson-reggeon vertex, we use the covariant representation [JJ Q21 [22]- For the 
production of two pseudoscalar particles (let it be im in the considered case), it reads: 



16vr2J+l A 271 " 1 

o = — - — -, olj = n — • 14 

The angular momentum operators are constructed of momenta p 1 - and k x which are orthogonal 
to the momentum of the two-pion system P = k\ + fc 2 : 



P P 1 1 



The coefficient £j normalizes the angular momentum operators, so that the unitarity condition 
appears in a simple form (for details see Appendix A). 

2.2.2 The t-channel 7r 2 -exchange 

The i?(7Tj)-exchanges dominate the spin flip amplitudes, and the amplitudes with m — 1 are 
here suppressed, see (JH]). However, their contributions are visible in the differential cross sections 
and should be taken into account. The effects appear owing to the interference in the two-meson 
production amplitude because of the reggeized ir 2 exchange in the t-channel. The corresponding 
amplitude is written as: 

(<*)+ 

5> Q/3 (W 2 ) -> nn)e^K 2 (s nN ,q 2 )^X%(k^) 9$\t) , (16) 

where ^7ri?(7r 2 ) — > 7T7t^ is the meson block of the amplitude related to the 7r 2 -reggeized 
t-channel transition, g^f 1 is the reggeon-pn vertex, H 7r2 (s 7rJ v, g 2 ) is the reggeon propagator, and 

e^g is the polarization tensor for the 2~ + state. Let us remind that k$ is the momentum of the 
outgoing nucleon. 

*i? = ^ »i? = 9„* - qj f- ■ (17) 

The 7r 2 particles are located on the pion trajectories and are described by a similar reggeized 
propagator. But in the meson block, the 2~ + state exchange leads to vertices different from 
those in the 0~ + -exchange, so it is convenient to single out these contributions. Therefore, we 
use for R^s^n, Q 2 ) the propagator given by ([TT]) but with eliminated 7r(0~ + )-contribution: 



Rn 2 (s, N , q 2 ) = exp (-z|4 Wm9) (9 2 )) 



(leadinq) / o\ 

\S-kNI S-kNO ) 



sin (f a% eadin9 \qZ)) T (\a { ! eadmg \q 



,2 
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Taking into account that 
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V e (a M a)+ - - ( a ±q a ±q + a ±g a ±q - -a ±g a ±q \ 

Z^ t af3 t a'f3' — 2 yyaa't/pp' <^ ypa'yaP' ^a^*'?' J ' 



(19) 



one obtains: 

-(2) f U U 



^a'(3'(^3 9 ) / ±q ±q , ±q ±q _ 2 ±q ±q \ _ 3 ^3a^3/f _ ^ m N ~ Q 2 ( _ a a a P 1 

9c 2 \9aa'9(3f3> + 9[3a'9af3' n9 a/3 9 a > j3> ~ 2 o „2 \ 9a/3 2 • 



In the large momentum limit of the initial pion, the second term in (J20l) is always small and 
can be neglected, while the convolution of k^ g k^ with the momenta of the meson block results 
in the term ~ s\ N . Hence, the amplitude for ^-exchange can be rewritten as follows: 



r, j^-Lq j^-Lq 

4~n angC) = -A a ^R^ 2 ) - t^JELM-R^^, q 2 ) (^Wft) 9^ . (21) 

A resonance with spin J and fixed parity can be produced owing to the ^-exchange with three 
angular momenta L = J — 2, L = J and L = J + 2, so we have: 

AtfMfa) - tttt) = 2 ^S(s)4S. w (p 1 )(- 1 ) J °tS , ( i 

J 

+E4 J) (^3i ^^Hltel 1 p)^U(^)^ 

J 

+ j2 p)x^u( fc± )6 . (22) 

The sum of the two terms presented in ( fit)]) and ( l2Ti) gives us an amplitude with a full set of 
the 7Tj-meson exchanges. 

Let us emphasize an important point: in the K- matrix representation the amplitudes 
A^.^s) (Eq. tfHD, j = leading, daughter-l) and A$(s), A { J \s), A { ^(s) (Eq. fl22])) 

differ only due to the prompt-production i^-matrix block (the term ^_R(t)(s) in (CQ)) while the 
final state interaction factor ([1 — p(s)K(s)]~ 1 in (JT])) is the same for each J. 

2.3 Amplitudes with aj-trajectory exchanges 

Here we present formulae for for leading and daughter ai-trajectories and leading (^-trajectory. 
2.3.1 Amplitudes with a^-trajectory exchanges 

The amplitude with t-channel ai-exchanges is a sum of leading and daughter trajectories: 

4^ CCt ° riCs) =T,A(irR(a?) - nn) (s nN , q 2 )i g^\t) , (23) 

«#> 
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where g^n 3 ^ is the reggeon-NN coupling and the reggeon propagator R a u) (s^jv, q 2 ) has the form: 



{SnN/SirNO) 



R 



i exp 



7T 

2' 



cos(faii ) (g 2 ))r(|ai J 1 ) (g 2 ) + |) ' 



(24) 



Recall that the a\ trajectories have a negative signature, £ n = — 1. Here we take into account 
the leading and first daughter trajectories which are linear and have a universal slope parameter 

ma Da GQi: 



pleading) ^ 



-0.10 + 0.72g 2 , a {daughte.r-l)^ _ _ 11Q + Qjty 2 . (25) 

As previously, the normalization parameter s n No is of the order of 2-20 GeV 2 , and the Gamma- 
functions in the reggeon propagators are introduced in order to eliminate the poles at q 2 < 0. 

For the nucleon-reggeon vertex we use two-component spinors in the infinite momentum 
frame, (p p and cp n , so the vertex reads (if^i(an z )ip p ) where n z is the unit vector directed 
along the nucleon momentum in the cm. frame of colliding particles. 

At fixed partial wave J PC = J ++ , the irR(a\) channel (j = leading, daughter-1) is charac- 
terized by two angular momenta L = J+1, L = J— 1, therefore we have two amplitudes for 
each J: 



A[iiR(a\ 



TXTX 



E 



,(-) 



.4 



(J+) 



(26) 



where the polarisation vector e 
the colliding particles as follows: 



,(-). 



(3 ~ n f3 , 



the GLF-vectors [23] defined in the cm. system of 



n 



(-) 



;i,o,o,-i)/2p 2 



;i,o,o,i)/2p 2 



(27) 



with p z — > oo. 



The products of Z and X operators can be expressed through vectors Vp and V^ J . 



j+i 



Pp 



-Pi 



y(J-) 



1 

J 



Pp 



-pi 



(28) 



Here, hi, pi and z are defined as: /c 2 = (k^k- 1 ), pi = (p L p L ),z = (—(k^-p^)^ / (^\f~ ^i\J ~ PlJ • 
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2.3.2 The amplitude with a 2 -trajectory exchange 

The amplitude with t-channel a 2 -trajectory exchange reads: 

£ (a)+ 

a S nN 

X%(k^)(^(k 3 )ij(p 2 )y^( q 2 ) , (29) 
where gffi is the reggeon-NN coupling and the reggeon propagator R a2 (s n N, <f) has the form: 

H-(^,, ) =exp )j sm(fttaa(92))r(X(92)) - PO) 

Recall that the leading a 2 trajectory has a positive signature, ^ = +1, it is linear with the 
following slope parameter [HI [19], [20] : 



a aa ( q 2 ) = 0.45 ± 0.05 + (0.72 ± 0.05)g 2 . (31) 

As previously, the normalization parameter s^no is of the order of 2-20 GeV 2 , and the Gamma- 
function in the reggeon propagator is introduced in order to eliminate the poles at q 2 < 0. 

Using Eqs. ([19"]). ([20]) . we obtain: 

A^ n = A a p(vR{a 2 )^)^\^R a2 {s« N ,q 2 ) U{k^{p 2 ) )g™tf) . (32) 

Due to a 2 exchange, the resonance with spin J can be produced from orbital momentum 
either J — 1 or J + 1. Thus, 



A a ,(nR(a 2 ) - tttt) = £ (A^( a )T<^ + ^( a )T^) , (33) 

j 



where 



tT ] = (34) 

Taking into account that the tensors T^p^ convolute with symmetrical tensor h^kf*, we 
obtain: 



0/3 30 3/3 " ^ j(j-i) r j(z \/=iH p ^ a/ — 



rp{J+l) h Lqi±q _ a J+l (Vgijl) ip// ,,±q,±q P±(J ~l) a J rr(J-l) u ±q u ±q foc\ 

K 3a K 3/3 — j £ pakP ^= Pf3 ^JK 2 ) K 3a K 3f3 (J+lWj x K 3a K 3f3 • \ 60 ) 
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2.3.3 Calculations in the Godfrey— Jackson system 

In the cm. system of the produced mesons, which is used for the calculation of the meson 
block (the GJ system), we write: 

(36) 

In this system the momenta follows: 

l£ p =p± = {0,0Ap), p 2 = (S + r f s ~ t)2 -ml, k*= S --ml, (37) 

k^ p = k± = (0, k sin cos (p, k sin sin if, k cos 0) , 
q = (g , 0,0, p), q = (s - ml + t) / (2y/s) , 

k 3 = (k 30 ,k 3x ,0,k 3z ), k 30 = (s nN - s - m 2 n )/(2^/s), k 3z = (2k 30 q - t)/(2p) . 



Recall that we use the notation A = (Ao,A x ,A y ,A z ) and cos© = z— —(k- L p ± )/(^J—k'^^J—p']_). 

4 J+) (hp-qp/2),V^ 



For the ai-exchange the convolutions Vg J+ ^ (k 3/3 — qp/2), vi 3 ' (k 3 p — qp/2) give us the am- 



plitude for the transition nR(ai ) into two pions (in a GJ-system the momentum k 3 is usually 
situated in the (x^)-plane). We write the amplitude in the form 

A^af) ^ ^ = J2 a JP J ~ lkJ (^ J \s)Y^e,^ + wl J \s)ReY}(e,^ , 

where the coefficients Wq J \s), W[ (s) are easily calculated: 

= E -Njo(h z -f] (\pfA^ -A^l ), 
'—^ y 2 J \ 7ra i * 7r7r 7ra i ^ 7r7r / 

< } = E-jn^hJ\p\ 2 JA^ +(J+1)A (J "! ). (38) 

i J(J+1) \ na\'-Hnr TTttj >7T7V J 



For a2-exchange, one has: 
T^k^ = ^P^h^-^NuImYKe,^) - ^k^fe,^)] (39) 



For the amplitude with orbital momentum J + 1, we write: 

TS +1) « = -W +1 ^(^-|)^/mF)(e,^) - ^^ y 1 ^^. (40) 
The final expression for the c^-exchange amplitude can be written as follows: 

Oh 

A&t^n = ^ E P J ~ lkJ [Wl 2 J Im Y) (O, + WZlmYftQ, <p)] x 

ZS 7rV J 

^a 2 (^,g 2 )(^(A;3)^(P2))^ 2) (g 2 ) , (41) 
13 



where 



W lJ 

a 2 



Nij fl 

—(to. - - T 



w, 



2.1 



N 2J k 
J 2 



3x 



«£4 ,(J) P "J 
J-l - 1 J+l +1 



(42) 



For the unpolarized cross section, the amplitude related to a 2 exchange does not interfere 
with either ir, ir 2 or a\ exchange amplitudes. If the highest moments are small in the cross 
section, one can assume that the combination in front of Y? is close to 0. Then 



W. 



2 



-N u (k 3z - ^) P 2 a J+1 A 



(J) 
+i ' 



W. 



2J 



0. 



(43) 



and, as a result, we have: 



np — >7rnn 



!¥- E yP J+1 k J [Nu(h z - P -)a J+1 A^ImY}(e, ip)] 

ZS ttN J J Z 



X R a2 (SnN, ?') ( ^(WM l^V) • 



(44) 



2.3.4 Partial wave decomposition 

The partial wave amplitude irR(a[^) — > tttt with fixed J ++ is presented in the i^-matrix form: 



a (l=j±i,j++) ( x _ k (l=j±i,j++). 



where K^ L J { f ) 1,J \s,q 2 ) is the following vector (b = tttt, KK, 7777, 7777', n-nix-n) 



■KR(aY'),b 



I -ip(s)K( J++ )(s)_ 



b.TVTT 



Here G 



K (L=J±l,J++) ( 2^ 



(L=J±l,J++,a)/ 2 



E 



G 



(L=J±l,J++,a) / 2 \ 



7rR(a< J,) ) 



(q 2 K J 



{J L=J±l, ++) 



Ml-s 
1 GeV 2 + srq \ s — sa 



s + s R o J s + sao 



(45) 



(46) 



(q 2 ) and & ^(g 2 ) are the g 2 -dependent reggeon form factors. 



2.4 7r j9 —> KKn reaction with KK- exchange by p-meson trajectories 

In the case of the production of a KK system the resonance in this channel can have isospins 
1 = and 1 = 1, with even spin (production of states of the types <fi and ao). Such processes 
are described by p-exchanges. 



11 



2.4.1 Amplitude with exchanges by p-meson trajectories 

The amplitude with t-channel p-meson exchanges is written as follows: 

47^r ies) = E4^(p;) "> KK ) R P ^N,q 2 )g^\ (47) 

where the reggeon propagator R Pj (s n N,q 2 ) and the reggeon-nucleon vertex gfn read, respec- 
tively: 

f i WJHi 2 ) 

D / 2n / - 71 " (j)f 2\\ [SirN S nN0 ) 

R fl {^ ) = e*p(-,-a«fe )j sin(fa Cfl (?2))r( i a « (?2) + l) . 

The pj-reggeons have positive signatures, £ p = +1, being determined by linear trajectories 
HEl CHI [20]: 

pleading) ^ ^ 0.50 + 0.83g 2 , a (f»«^-i) ( q 2 ) ~ -0.75 + 0.83g 2 . (49) 

The slope parameters are in (GeV/c)~ 2 units, s^ato ~ 2 — 20 GeV 2 . Two vertices in ^ cor- 
respond to charge- and magnetic-type interactions (they are written in the infinite momentum 
frame of the colliding particles). 

The meson-reggeon amplitude can be written as 

A^Rfa) - KKj = E e *<->1*W^(PM > ( 5 °) 

where the polarisation vector eL was introduced in fl36l) . 

We use the convolution of the Z and X operators in the GJ-system (see notations in ( 1371) : 



j 

The convolution of the spin-momentum operators in ( 1501) gives: 



.(51) 



A ( m - vttt) = £ Zfp'k'htoy/iNjJm Y}(Q, (p)A^ t) KR (s). (52) 
Let us remind that in the GJ-system the vector fc 3 is situated in the (x^)-plane. 
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2.4.2 Partial wave decomposition 



The amplitude for the transition irR Pj (q 2 ) — > KK in the i^-matrix representation reads: 



.4 



where K^ R ^ b (s, q 2 ) is the following vector {b = irn, KK, rjrj, rjrj' ', Trnnn): 



(J++) 

irR(pj),KK 



» = Y. K i J R{ P ])A S ^ 



I -ip(s)K( J++ )(s)\ b 



:KK 



K (J++) (s a 2 ) 



G 



(J++,a) f 2 \ (J ++ ,a) 
■R( Pj ) 



Wb 



+ F 



a 1 y± a a 

(J ++) 2 1 GeV 2 + sub's s - s A 



TrR{pj),b 



(Q 2 



S + S R0 



- 



s + sao 



Here G^ R ^\q 2 ) and F^ R ^ b (q 2 ) are the reggeon g 2 -dependent form factors. 



(53) 



(54) 



3 Low-energy three-meson production in the i^-matrix 
approach 

Here we present elements of the i^-matrix technique for the low-energy reactions pp — > mrir, Trrjr), 
7rKK. The .ff -matrix technique provides a compact and, hence, a convenient way for studying 
resonances in multiparticle processes of such a type. However, we have to pay a price for the 
simplifications the iT-matrix technique gives us: we cannot take into account in a full scale 
the partial wave left singularities as well as the singularities related to the rescattering of all 
particles (for example, the singularities of the triangle type diagrams) 

The use of the i^-matrix approach to the combined analysis of different processes is based 
on the fact that the denominator of the i^-matrix two-particle amplitude, [1 — f>K]~ l is common 
for all processes, depending only on quantum numbers of the considered two-meson system. 

Let us illustrate this statement using as an example the amplitude of the pp annihilation from 
the 1 Sq level: pp^So) — > threemesons. In the .fT-matrix approach, the production amplitude 
for the resonance with the spin J = in the channel (1 + 2) reads: 

A 3 (s 12 ) ca = Y,(4 pr ° mpt \si2)) ( - - L . . ) , (55) 

b V / cb\\ - ipi2Ki2{S 12 )S ba 

where c = pp^So) and a,b E inr,r)r), KK. The denominator [1 — ipviKi2{si'i)]~ 1 depends on 
the invariant energy squared of mesons 1 and 2 and it coincides with the denominator of the 
two-particle scattering amplitude. The factor K^ )Tompt \si2) stands for the prompt production 
of particles and resonances in this channel: 

-(prompt) f \\ ST^ c iJb 



Kr° mpt, (s 12 ) =£^- + ^12), (56) 

'cb n K~ s 12 
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where A( n ) and ip c b are the parameters of the prompt-production amplitude, and and /x n 
are the same as in the two-meson scattering amplitude. 

The whole amplitude for the production of the (J = 0)-resonances is defined by the sum of 
contributions from all channels: 

A 3 (s 12 ) + A 2 (s 13 ) +A 1 (s 23 ). (57) 

The amplitudes A 2 (s 13 ) and A^s^) are given by formulae similar to fl55|) . fl56l) but with different 
sets of final and intermediate states. 

To take into account the resonances with non-zero spins J, one has to substitute in fl55|) 




where the i^-matrix amplitude A 3 (si 2 ) is determined by an expression similar to (155]) . 

The analysis performed in [Ml 125] showed that in the reactions pp(ai rest) — > 7r 7r 7r°, ir°7r°r), 
7r°?7?7 the determination of parameters of resonances produced in the two-meson channels does 
not require the explicit consideration of the triangle diagram singularities — it is important 
to take into account only the complexity of parameters Aj™ and tp a b in ( |56l) which are due 
to multiparticle final-state interactions. Note that this is not a universal rule for the meson 
production processes in the pp annihilation - for example, in the reaction pp — > 7]7t + 7t~7t + tt~ 
|26j . the triangle singularity contribution is important. 

4 Fitting procedure 

Here we present results of the combined fit for low energy annihilation reactions pp — > tttttt, 
mrr], Trrjr) and high energy peripheral production ir~p — > tt°it + n. 

4.1 The K- matrix fit of annihilation reactions at rest pp into 7T7r7r, 

7T7TT], TTTJIJ 

We have included into the fit procedure the following data sets for the production of three 
mesons in pp annihilation: 

(1) Crystal Barrel data onpp(atrest, from liquid H 2 ) — > 7r°7r 7r°, 7r°7r°?7, 7r°r]r] [27] and 

(2) the data in gas pp(atrest, from gaseous H 2 ) — > 7r°7r 7r°, 7r°7r°?7 [28] [29]. 

The considered i^-matrix amplitude takes into account innnr, KK and 777/ channels as well 
- parameters for these channels are taken from [TO] . 

First, we present the formulae for the reactions pp — > 7t°it 7t , 7r°7r°?7, 7r r/?7 from the liquid 
H 2 , when annihilation occurs from the 1 S$pp state and scalar resonances, fo and do, are formed 
in the final state. This is a case which represents well the applied technique of the three-meson 
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production reactions. A full set of amplitude terms taken into account in the analysis [TU] 
(production of vector and tensor resonances, pp annihilation from the P-wave states 3 Pi, 3 P 2 , 
1 Pi) is constructed in an analogous way. 

(i) Production of the S-wave resonances. 

For the transition pp ( 1 S'o) — > 7r 7r°7r° with the production of two pions in a (00 ++ )-state, 
we use the following amplitude: 



A 



pp (l 1 So)— >7r 7r°7r 



(59) 



App (l 1 S o )7r0,7r07r o ( s 23) + A pp (V- S )ir° ( s 13) + A pp (li5 )7r ,7r 7r° ( s 12) 



The four-spinors VK — P2) and ijj(pi) refer to the initial antiproton and proton in the I^ 2S+1 ^Lj = 
l 1 ^ state. For the produced pseudoscalars we denote amplitudes in the left-hand side of (|59p 
as A pP (ii So )p e (sij)- 

The amplitudes for the transitions pp (O 1 ^) — > r77r°7r°, pp (l 1 ^) — * ir r]r] have a similar 
form: 



and 



App (OiS )-*riTT -K 



App (Ql Sq)t],tt 7t { s 23) + App (0 1 So)7r°,^7r°( s 13) + A pp (p 1 S )-ir ,rr*° ( S 12, 



(60) 



App (iis )^ir°rm - ^(-^2); 



«75 



-1p{p 1 



App (ll5 )7r ,r ) 7 ? ('523) + ^4pp (l 1 ^)^,^ ( s 13) + A pp (±1 S )r],rin (Sl2) 



(61) 



For the description of the S'-wave interaction of two mesons in the scalar-isoscalar state (index 
(00)) the following amplitudes are used in ( l59l) . ( l60l) and ( JBTT) : 



ab 



(62) 



Here 6 = 7r°7r°, 7777 and a = 7r°7r°, 7777, i^^, 7777', 7r°7r 7r 7r°. The i^-matrix term is responsible 
for the two-meson scattering. The i^-matrix terms which describe the prompt resonance and 
background meson production in the pp annihilation read: 



- fl pp(l 1 S , )T°,al 6 23j 



E 



A(00,a) (a ) 
1V pp(l 1 S )n°ya 



Ml - S 



23 



+ 



(oo) 1 GeV + s \ ( s 23 - s A 

^(liSoW 



S23 + s 



S23 + sao 



(63) 
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The parameters A^'^^o a and ^p^ig^o a are complex-valued, with different phases due to 
three-particle interactions. Let us recall: the matter is that in the final state interaction term 
we take into account the leading (pole) singularities only. The next-to-leading singularities are 
accounted for effectively, by considering the vertices pp — > mesons as complex-valued factors. 

(ii) Three-meson amplitudes with the production of spin-non-zero resonances. 

In the three-meson production processes, the final-state two-meson interactions in other 
states are taken into account in a way similar to what was considered above. 

The invariant part of the production amplitude A^ty 1Sq fe )(23) for the transition (/ 1 5'o) — > 
1 + (2 + 3) t j, where the indices tj refer to the isospin and spin of the meson in the channel 
b = 2 + 3, is as follows: 



a 

a(*7» (a) 



all 



s 23 



The parameters ^^'nis )v ^wj 1 s )i a ma y be complex- valued, with different phases due to 
three-particle interactions. 

The .fT-matrix elements for the scattering amplitudes (which enter the denominator of (l64"j) ) 
are determined in the partial waves 02 ++ , 10 ++ , 12 ++ as follows: 

(1) Isoscalar-tensor, 02 ++ , partial wave. 

The -D-wave interaction in the isoscalar sector is parametrized by the 4x4 i^-matrix where 
1 = 7T7T, 2 = KK, 3 = nn and 4 = multi — meson states: 

I<T\s) = DM |Jg + /r 1G ;^ S j D b (s) . (65) 

Factor D a (s) stands for the .D-wave centrifugal barrier. We take this factor in the following 
form: 

D 'W = kT^ a > a=1 ' 2 ' 3 ' ( 66 ) 



where k a = Js/A — m\ is the momentum of the decaying meson in the cm. frame of the 
resonance. For the multi- meson decay the factor D^(s) is taken to be 1. The phase space 
factors we use are the same as those for the isoscalar S'-wave channel. 

(2) Isovector-scalar, 10 ++ , and isovector-tensor, 12 ++ , partial waves. 
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For the amplitude in the isovector-scalar and isovector-tensor channels we use the 4x4 
if-matrix with 1 = Tirj, 2 = KK, 3 = irt]' and 4 = multi- meson states: 

AfiOW = D.(s) + jfiflH^!±fi) 0t(s) . ( 67 ) 

Here j = 0,2; the factors D a (s) are equal to 1 for the 10 ++ amplitude, while for the .D-wave 
partial amplitude the factor D a (s) is taken in the form 



2 



D a (s)= U2 ^%, 2 , a = 1,2, 3, D 4 (s) = l. (68) 
The results of the fit for the two solutions discussed below are shown in Figs. 



4.2 The A'-matrix fit of high-energy meson production: the tt- and 
ai-trajectory exchanges 

The leading contribution from the 7r-exchange trajectory can contribute only to the moments 
with m — 0, while the ai-exchange can contribute to the moments up to m = 2. The char- 
acteristic feature of the ai exchange is that moments with m = 2 are suppressed compared to 
moments with m = 1 by the ratio k^ x /k^ z which is small for the system of two final mesons 
propagating with a large momentum in the beam direction. 

The amplitudes defined by the tt and a± exchanges are orthogonal if the nucleon polarisation 
is not measured. This is due to the fact that the pion trajectory states are defined by the 
singlet combination of the nucleon spins while the a\ trajectory states are defined by the triplet 
combination. This effect is not taken into account for the S-wave contribution in (Q which can 
lead to a misidentification of this wave at large momenta transferred. 

The 7r 2 particle is situated on the pion trajectory and therefore should be described by the 
reggeized pion exchange. However, the 7r 2 -exchange has next-to-leading order contributions 
with spherical functions at m > 1. The interference of such amplitudes with the pion exchange 
can be important (especially at small t) and is taken into account in the present analysis. 



4.3 Results of the fit 

To reconstruct the total cross section of the reaction n~p — > 7r°7r n which is not available to us 
we have used two partial wave decompositions provided by the E852 collaboration [TTJ . The 
cross section was reconstructed by Eq. (Q and decomposed over moments. The two partial 
wave decompositions produced very close results for the moments and we included the small 
differences between them as systematical errors. 

The 7i~p — > 7i°7T°n moments can be described successfully with only n, a% and 7r 2 leading 
trajectories taken into account and a simple assumption about the t-dependence of form factor 



20 



for all partial waves. Moreover, we have found two solutions which differ by their contributions 
from these exchanges. Such an ambiguity is likely to be connected with the lack of polarisation 
data and can be resolved by data from future experiments. 

The quality of the description of the Crystal Barrel data by both solutions is shown in Figs. 
I2"iri"3l and the x 2 is given in Table [TJ Here we also provide x 2 values for the tttt — > rjr] and 
7i7i — > rjr]' S-wave extracted by the GAMS collaboration [12], [13] from the tt~N data taken at 
small transfer energies. It is seen that both solutions describe Crystal Barrel and GAMS data 
with the same quality. The main reason that the K-matrix parameters for the S and D- waves as 
well as P- vectors for the pp annihilation into three mesons are very similar in the two solutions. 
However the P-vectors for description of the E852 data are different in solutions 1 and 2. 

In a more detail: these two solutions differ by the fraction of the tt, a\ and ir 2 exchanges 
already in the region of small energy transferred. The first solution has a very large, practically 
dominant contribution from the a\ exchange to the D-wave (see Fig. [TJJ) . The contribution 
from the a\ exchange to the S'-wave is small. In this solution there is no notable signal from 
the /q(1300) state either at small or at large energy transferred. If /o(1300) is excluded from 
this solution, only the description of the Crystal Barrel and GAMS data is deteriorated while 
the description of the E852 data has the same quality. 

In the second solution the D-wave at small energies transferred is dominantly produced 
from the tt exchange. The fraction a\ exchange at \t\ < 0.1 is about 2.5 — 3%. At large 
energy transferred, like in solution 1, the contribution from a% exchange becomes comparable 
and even dominant. The S-wave has a well known structure at small \t\. At intermediate 
energies the contribution from the a\ exchange becomes dominant and a signal from the /o(1300) 
state is well seen in this wave. At very large \t\ (—1.5 < t < —0.4 GeV 2 ) the contribution 
from a\ exchange is rather small. The dominant contribution comes from the /o(980) state 
produced from tt exchange. Here our analysis is in contradiction with the result reported by the 
E852 collaboration which observed a strong S-wave signal around 1300 MeV in this t-interval. 
However, the contribution from / (1300) at intermediate energies transferred is important for 
the description of data with this solution. If this state is excluded from the fit the description 
is notably deteriorated, see Tabled], solution 2(-). This subject is considered in the following 
section in detail. 

The Krakow group reported from the analysis of the polarized data that at small t the 
dominant contribution comes from the 7r-exchange [30]. They point out that the second solution 
is possibly a physical one. However, the final conclusion can be made only after including these 
(yet unavailable to us) data in the present combined analysis which uses reggeon exchanges. 

The description of the moments at small and large \t\ for the two solutions is shown in 
Figs. [TB] and [171 correspondingly. The second solution produces a systematically better overall 
description except for the Y± moments at large energies transferred. 

The S-wave was fitted to 5 poles in the 5-channel .fT-matrix, described in detail in previous 
sections. The parameters for the first solution are very close to those for the second one, e.g. 
the parametrization given in Table [2] describes both solutions, and the given errors cover a 
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Table 1: The x 2 per data point for the description of the Crystal Barrel and GAMS data. Two 
solutions are given as well as that with / (1300) excluded from the fit, solution 2(-). For the 
E852 data the x 2 is calculated for all moments in the given t-interval 



Data 


Solution 1 


Solution 2 


Solution 2(-) (no / (1300)) 


PP ~ > 


7r°7r°7r (Liq) 


1.360 


1.356 


1.443 


pp — ► 


7T 7r°7r (Gas) 


1.238 


1.242 


1.496 


pp — > 


7]7r°r] (Liq) 


1.350 


1.442 


1.446 


pp — > 


?77r ?7 (Gas) 


1.503 


1.371 


1.315 


pp — > 


ttV° (Liq) 


1.210 


1.236 


1.412 


pp — > 


7rV°(Gas) 


1.099 


1.119 


1.227 


nil — 


• 7777 (S-wave) 


1.08 


1.19 


1.38 


7T7T — 


• 7777' (S-wave) 


0.26 


0.41 


0.45 



marginal change in both descriptions. 

The D-wave was fitted to 4 poles in the 5-channel (inr, KK, 7777, uou and Att) i^-matrix. The 
position of the first two D-wave poles was found to be 1270 — 797 MeV and 1530 — i72 MeV 
which corresponds to the well-known resonances ^(1270) and / 2 (1525). The third state has 
a Flatte-structure near the ouu threshold and is defined by two poles on the sheets defined by 
the ujoj cut. Due to the fact that we do not fit directly the ujoj production data these positions 
can not be defined unambiguously. For example, in the framework of the solution 1 (dominant 
ai-exchange in the D-wave) we found at least two solutions for the pole structure in the region 
of 1560 MeV. In the first the pole is situated at 1565 — 7140 MeV on the sheet above the ujuj 
threshold and 1690 — 7290 MeV on the sheet below the ujuo threshold. In the other solution 
the position of the pole is 1530 — 7262 and 1699 — 7216, correspondingly. The closest physical 
region is for both poles the beginning of the ujuj threshold M ~1570 MeV, where they form 
a relatively narrow (220-250 MeV) structure which is called the /2(1560) state, see Fig. [TSJ 
A similar situation was observed in the solution 2. The K-matrix D-wave parameters for the 
solution 2 are given in Table [31 

The fourth D-wave i^-matrix pole, /| are (1980) cannot be rigidly fixed by the present data. 
The position of the corresponding pole is also not stable: one can easily increase the mass of 
the pole with the simultaneous increase of the width, spoiling only slightly the description of 
data. Because of that we consider this pole as some effective contribution of resonances located 
above 1900 MeV. 

The 7T7T — > 7T7T S-wave elastic amplitude for the second solution is shown in Figfl9l The 
structure of the amplitude is well known, it is defined by the destructive interference of the 
broad component with / (980) and / (1500). Neither / (1300) nor / (1750) provide a strong 
change of the amplitudes. However, this is hardly a surprise: both these states are relatively 
broad and dominantly inelastic. 

The 7T7T — > 7T7T D-wave elastic amplitude is shown in Fig. [20j The amplitude squared is 
dominated by the ^(1270) state. Neither of /2(1560) and ^(1510) (which are included into 
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Table 2: Masses and couplings (in GeV units) for the S-wave .fT-matrix poles (/o are states) as 
well as the amplitude pole positions (given in MeV). The II sheet is defined under the mr and 
47r cuts, the IV sheet is under hit, 4ir, KK and 7777 cuts, and the V sheet is determined by tttt, 
4ir, KK, rjr) and 777/ cuts. 





a = 1 


a = 2 


a = 3 


a = 4 


a = 5 


M 


0.7201;™, 


-1 r,nn+- 040 
l.ZZU„ 030 


1.210 ±030 


J..OOU_ 020 


1.850 ± .040 


(a) 
% 


0.7601;°™ 


0.820 ± 0.090 


0.470 ± .050 


0.360 ± .050 


0.440 ± .050 


(a) 

9h 








0.850 ± .100 


0.570 ± .070 


-0.900 ± .070 




-(60 ±12) 


28 ± 12 


30 ± 14 


8± 15 


-(52 ± 14) 




a = TTTT 


a = KK 


a = 


a = 7777' 


a = 47T 


fla 


0.180 ± .120 


0.150 ± .100 

/6a = 


0.240 ± .100 

b = 2,3,4,5 


0.300 ± .100 


0.0 ± .060 


II 

sheet 


1030 + fo 
-*(35!Jg) 




Pole position 






III 

sheet 


850^ 
-i(100±25) 


IV 

sheet 




1290 ± 50 
-i(170^g) 


1486 ± 10 
-i(57±5) 


1510 ± 130 

-z(800l^) 




V 

sheet 


1800 ± 60 
-«(200 ± 30) 
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Table 3: Masses and couplings (in GeV units) for D-wave i^-matrix poles (/| are states) for the 
solution 2. The III sheet is defined by irn and 4ir and KK cuts, the IV sheet by tttt, Att, KK 
and ujuj cuts. The values marked by * were fixed in the fit. 





a = 1 


a = 2 


a = 3 


a = 4 


M 


1.286 ±0.025 


1.540 ±0.015 


1.560 ±0.020 


Z.ZUU_g 200 


J7T7V 


0.920 ±0.020 


-0.05 ±0.080 


0.280 ±0.100 


-0.30 ±0.15 


g(«) 
i/rm 


0.420 ± 0.060 


0.27 ±0.15 


0.400 ±0.200 


1.2 ±0.6* 


(2) 


-0.150 ±0.200 


0.370 ±0.150 


1.170 ±0.450 


1.0 ±0.4 


o(«) 


0* 


0* 


0.540 ±0.150 


-0.05 ± 0.2 




a = TTTT 


a = f]7] 


a = ujuj 


a = 4ir 


fla 


0.03 ±0.15 


-0.11 ±0.10 


0* 


0* 


f2a 


-0.11 ±0.10 


-1.8 ±0.60 

fba = 


0* 

6 = 3,4,5 


0* 






Pole position 




III sheet 


1.270 ±0.008 
-i 0.097 ±0.008 


1.530 ±0.012 
-i 0.064 ±0.010 






III sheet 






1.690 ±0.015 
-i 0.290 ±0.020 




IV sheet 






1.560 ±0.015 
-i 0.140 ±0.020 





the KK channel of the i^-matrix) show a meaningful structure in the amplitude squared. The 
i^-matrix parameters found in the solution are given in Table [3j 

4.4 The / (1300) state 

In the solution 2 the fit of the E852 data shows a large contribution from the /o(1300) state to 
Yq moment due to ai exchange at —0.2 < t < —0.1 and —0.4 < t < —0.2 GeV 2 . At very small 
(—0.1 <t< —0.01 GeV 2 ) and large (— 1.5 <t<— 0.4 GeV 2 ) energy transferred the contribution 
of this state to the Y" ° moment is less pronounced. If the K-matrix pole which corresponds to 
the /o(1300) state is excluded from the fit (all couplings are put to zero) the total \ 2 changes 
rather appreciably. The corresponding values for the description of the Crystal Barrel and 
GAMS data are given in the last column of Tabled] (solution 2(-)). The mass slices made in 
the region of the state show systematical discrepancies in this case (see Figj2"T]). 

For the description of the E852 data the main effect is seen, as expected, for the second 
and third t intervals. The comparison of the solutions with and without /o(1300) for these t 
regions is shown in Fig. [521 Here the description of the Y± moment is systematically worse for 
the fit where /o(1300) is excluded. The x 2 P er data points change for this moment from 1.84 to 
3.63 for the -0.2<t<-0.1 GeV 2 interval and from 2.07 to 4.90 for the -0.4<t<-0.2 GeV 2 
interval. The fit without /o(1300) produces a worse description also for Y® and V 4 °. At intervals 
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of small and large t the description has the same quality and can hardly be distinguished on 
the pictures. The contribution of the S-wave to the moment Yq from this solution is shown 
in Fig. [23j It is seen that an appreciable contribution from a\ exchange at the mass region 
1300-1500 MeV is needed and the fit tries to simulate it (although not very successfully) by an 
interference between the broad component and the /o(1500) state. 

Below we present the pole positions of the S-wave amplitude (in MeV units) and couplings 
calculated as pole residues (in GeV units): A a ^ b ~ G a Gb[(M — iT/2) 2 — s]" 1 + smoothterms 
with a,b — 7T7T, KK, rjr], rjr]', irinnr; the couplings are written as G a = g a exp(iip a ), the phases 
are given in degrees. For resonances / (980), / (1300), /o(1500), / (1200 - 1600), / (1750) we 
obtain: 





/o(980)i stpo / e 


/o(980) 2n( ipo/e 


/o(1300) 


/o(1500) 


/o(1200 - 1600) 


/o(1750) 


M 


1030±f8 


85011° 


1290±50 


1486±10 


1510±130 


1800±60 


r/2 


oo_ 16 


100±25 


1701$ 


57±5 


8001}°° 


200±30 


Sheet 


II 


III 


IV 


IV 


IV 


V 




0.42±0.07 


0.39±0.05 


0.28±0.08 


0.24±0.05 


0.82±0.06 


0.55±0.05 




-71±8 


45±7 


27±10 


65±8 


10±12 


i0 -15 


9kk 


0.62±0.06 


0.68±0.12 


0.15±0.05 


0.17±0.04 


0.84±0.08 


0.11±0.04 




3±8 


155±6 


35±15 


48±8 


2±10 


55±20 


9r\Tj 


0.51±0.07 


0.58±0.10 


0.14±0.06 


0.10±0.03 


0.40±0.06 


0.18±0.05 




10±8 


157±10 


57±8 


96±6 


16±7 


40±12 




0.42±0.08 


0.46±0.12 


0.17±0.07 


0.18±0.06 


0.14±0.05 


0.35±0.07 




18±8 


160±10 


75 ±15 


143±15 


80±17 


18±6 


9iir 


0.16±0.05 


0.29±0.10 


0.80±0.15 


0.47±0.08 


1.30±0.20 


0.85±0.20 




25i? 5 


155±12 


205±12 


156±10 


5±12 


150±14 



(69) 



4.5 Isovector scalar and tensor resonances below 1.7 GeV 

The isovector states are contributing strongly to the pp — > rjrjTT and pp — > 7r°7r°?7 amplitude. 

The ao(980) is clearly seen on the 777771"° Dalitz plot and on the 7r77 mass projection. A success- 
ful description of these data can be obtained with the Flatte parametrization of this resonance 
(pole parametrization with decays into channels irrj and KK). Within this parametrization 
(for more detail see [21]) we obtain the following masses and couplings for a (980): 

M (MeV) <^(GeV) R = g 2 KR /gl v (7()] 
986 ±4 0.175 ±0.015 1.20 ±0.15 1 ; 

These parameters result in the following two poles on the I-st sheet (under the irrj cut) and the 
Il-nd sheet (under the irrj and KK cuts) for ao(980): 

I-st sheet Il-nd sheet . 

1000 ± 6 — z(35 ± 4) MeV, 940 ± 20 — i(85 ± 15) MeV ^ ' 



25 



The second isovector scalar state ao(1474) [T] (defined as ao(1450) in PDG) is situated in 
the 1500 MeV region. We found: 



Pole position of a ( 1474) : M = 1515 ± 30 - i(115 ± 18) MeV. (72) 

This state is highly inelastic with a branching ratio into the nr) channel less than 10%. We 
have not found any indications for an extra isovector scalar state in the mass region between 
a (980) and a (1474). 

The resonance a 2 (1320) contributes to the pp — *> 7r°7r°?7 strongly, here we see also a 2 (1675)- 
signal. We parametrized these resonances in the Breit-Wigner form and have found the following 
amplitude poles: 

a 2 (1320) a 2 (1675) . , 

M = 1309 ± 4 — z'(55 ± 2) MeV, M = (1675 ± 25) - i(135lfo) MeV 1 } 

For a2(1320) we have found a mass 1309 ± 4 MeV which is lower by 9 MeV than the average 
PDG value but corresponds very well to the analysis of high statistical data performed by the 
VES collaboration [39]. The observed width of the state 111 zb 4 corresponds well to other 
observations from the Tin decay mode. 

The a 2 (1675) state improves the fit; however, the mass and width of this state can not be 
well defined from these data because the resonance is situated on the phase volume boundary 
and is suppressed by the D-wave centrifugal barrier. The obtained values are compatible with 
previous findings of the Crystal Barrel collaboration [34]. Let us mention that the mass and 
width of this state can be much better defined from the L3 data on 77 interaction into 7r + 7r - 7r° 



5 Summary for isoscalar resonances 

We develop a method for the analysis of the reactions ttN — > two mesons + N at large energies 
of the initial pion. The approach is based on the use of the reggeized exchanges that allow 
us to analyze simultaneously the data obtained at small and large momentum transfers. In 
the present article the method is applied to the analysis of the tt~N — > N data measured 
by the E852 experiment. The inclusion of the Crystal Barrel data on the proton- antiproton 
annihilation at rest into the 37r°, 7r°?7?7 and 7r°7r°?7 channels helps to reduce ambiguities in the 
isoscalar sector and investigate the properties of the isovector scalar and tensor states. 

As the result of the analysis the K-matrix parameters of the isoscalar-scalar and isoscalar- 
tensor states was obtained up to the invariant mass 2 GeV and pole positions of corresponding 
amplitudes are defined. 
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5.1 Isoscalar-scalar sector 



In the scalar sector the contribution of the / (1300) is necessary to get a consistent description 
for the data set analyzed: 

Pole position of /„( 1300) : M = 1290 ± 50 - i(170+|°) MeV. (74) 

According to our fit, the strong signal in the 7T7T spectrum in the region 1300 MeV is formed 
by two contributions, by / (1300) (dominantly the a\ reggeized exchange) and / 2 (1275) (the 7r 
and a\ reggeized exchanges). 

The position of the / (980) is defined very well. The resonance reveals a double pole 
structure around the KK threshold. 

Pole positions of /„ (980) : (75) 
sheet II (under vrvr and ixixixix cuts) : M = 1030±fg - z(35±^) MeV, 
sheet III (under mi, nnixix and KK cuts) : M = 850l|] - «(200 ± 50) MeV. 

The /o(1500) is defined from the combined fit with a good accuracy: 

Pole position of /„( 1500) : M = 1486 ± 10 - z(57± 5) MeV. (76) 

The broad state /o(1200 — 1600) (the scalar glueball descendant) gives contribution in titi 
scattering amplitudes in region up to 2 GeV; the following pole position is found 

Pole position of /„( 1200 - 1600) : M = (1510 ± 130) - i(800+^°) MeV. (77) 

The / (1750) is a dominantly ss state [1] and is needed to describe rnr — > nn and nil — > 7777 
amplitudes above 1750 MeV. 

Pole position of / (1750) : M = 1800 ± 60 - i(200 ± 30) MeV. (78) 

Parameters of this state differ from that observed by the BES [36] and WA102 [37] collaborations 
(denoted as /o( 1710)); one should, however, have in mind that in the case of strong interferences 
characteristics of a peak in the data does not correspond to the resonance position. A combined 
fit of the Crystal Barrel, CERN-Munich, E852, GAMS and BES data is needed and it is one 
of our future objectives. 



5.2 Isoscalar-tensor sector 



The D-wave reveals the resonances ^(1275), ^(1525), /2(1565), and /2(1950) with the following 
pole positions: 



/ 2 (1275) 
/ 2 (1525) 
/ 2 (1565) (2nd solution) 



M = 1270 ± 8 - i(97 ± 8) MeV, 
M = 1530 ± 12 - z(64 ± 10) MeV, 
Mj = 1690 ± 15 - z(290 ± 20) MeV, 
M n = 1560 ± 15 - i(140 ± 20) MeV, 



(79) 
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In the case of /2(1565) the K-matrix fit can be obtained only with the large coupling of this 
state to ujuj (and, possibly, to pp) channel (note that this result is in a very good agreement 
with the analysis of the proton-antiproton annihilation into ujutt [38] ). The large coupling to 
ujuj leads to the double pole structure of /2(1565), see Fig. [18j 

The state / 2 (1980) can not be identified unambiguously from the present data due to its 
large inelasticity. It plays the role of some broad contribution needed for the description of the 
ttN data 

5.3 Isoscalar sector J PC = 4 ++ 

For the description of high moments in the ttN — > ir 7r°N data a contribution from a 4 ++ state 
is needed. This state is identified as /4(2025). Due to the lack of data at high masses this state 
was fitted as a two channel (7T7T and 4ir) one pole K-matrix. 



Here, as previously, masses and couplings are in GeV units. The position of the pole is equal to 
(1966 ± 25) — % (130 ± 20). The amplitude phase and the Argand diagram for the isoscalar 4 ++ 
state is shown in FigfSH The tttt — > nir 4 ++ amplitude has a peak at 1995 MeV and is slightly 
asymmetrical: the half height is reached at the mass 1880 and 2165 MeV. The branching ratio 
of the 7T7T channel at the pole position is 20 ± 3% which is in agreement with the PDG value 
within the error. 
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6 Appendix A. Angular Momentum Operators 

The angular-dependent part of the wave function of a composite state is described by operators 
constructed for the relative momenta of particles and the metric tensor. Such operators (we 
denote them as Xj£> „ , where L is the angular momentum) are called angular momentum 
operators; they correspond to irreducible representations of the Lorentz group. They satisfy 
the following properties: 

(i) Symmetry with respect to the permutation of any two indices: 



M (GeV) g n% g 4n Un^nn 

1.970 ±30 0.550 ±0.050 0.490 ± 0.080 -0.025 ± 0.050 



(80) 



fll...Hi...flj...flL 



XV 



(81) 
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(ii) Orthogonality to the total momentum of the system, P = k\ + k 2 : 

P^l,,.., L = 0. (82) 

(iii) Tracelessness with respect to the summation over any two indices: 

Qa u X^ „ „ „ = 0. (83) 

H^t*] fii...fj,i...fj,j ...fi L \ J 

Let us consider a one-loop diagram describing the decay of a composite system into two spinless 
particles, which propagate and then form again a composite system. The decay and formation 
processes are described by angular momentum operators. Owing to the quantum number 
conservation, this amplitude must vanish for initial and final states with different spins. The 
S-wave operator is a scalar and can be taken as a unit operator. The P-wave operator is a 
vector. In the dispersion relation approach it is sufficient that the imaginary part of the loop 
diagram, with S- and P-wave operators as vertices, equals 0. In the case of spinless particles, 
this requirement entails 

4^ = > ( 84 ) 

where the integral is taken over the solid angle of the relative momentum. In general, the result 
of such an integration is proportional to the total momentum P^ (the only external vector): 

Convoluting this expression with P^ and demanding A = 0, we obtain the orthogonality con- 
dition (1821) . The orthogonality between the D- and S- waves is provided by the tracelessness 
condition flHBl : equations (1821) . flHBl provide the orthogonality for all operators with different 
angular momenta. 

The orthogonality condition fl82l) is automatically fulfilled if the operators are constructed 
from the relative momenta and tensor g^ v . Both of them are orthogonal to the total 
momentum of the system: 



l \ / i , N I PiiPiJ 

T 



k i = ^9i„(h - k 2 )u , gi u = g^u - • (86) 



In the cm. system, where P = (Pq, P) = {y/s, 0), the vector k is space-like: k L = (0, k, 0). 

The operator for L = is a scalar (for example, a unit operator), and the operator for L = 1 
is a vector, which can be constructed from only. The orbital angular momentum operators 
for L = to 3 are: 

X(°)(^) = l, XjP = k£, (87) 



v-(3) r k ±\ _ 5 r, ± ,±,± _ k± / i , ± i , ± i , j_ 
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The operators Xj£> for L > 1 can be written in the form of a recurrency relation: 



L 



W-/*L\ / £2 Ml-Mi-lMt+l-MlA ' y Mi" 

2 L 



2L 



i iJ.i—iJ.i-iiJ.i+i—iJ-j-itMj+i—iJ-Lciv v ))■ 

1 i,,=l 



The convolution equality reads 

^J^^^^fciLM- (89) 

On the basis of Eq. (1891) and taking into account the tracelessness property of Xj£> , one can 
write down the orthogonality-normalisation condition for orbital angular operators 

A 2/ - 1 

a L = Y[—j— (9°) 
Iterating equation (|88|) . one obtains the following expression for the operator Xj£> flL : 



Mi—Mi 

2L- 1 



J.-L J.-L jU-L k- 1 - jU-L 

Ml M2 M3 M4 Mi 



%M2>3>4 ' ' ' Vi ^ "MlM3 fsV* ' ' ' Mi ^ • • • 



K 1 1.-L 

Mi 



J i r,- 1 - r,- 1 - t- 1 - t- 1 - t 

T ^2L — 1)(2L — 3) x^^^y^*, Ms M6 

M1M2 "M3 Ms M4 M6 • • • Vi 



(91) 



The projection operator O^;;;^ (or 0*£"££(A. P)) is constructed of the metric tensors 5^,. 
It has the properties as follows: 

v-(i) /QMi-Mi _ 
^ v Mi — Mi fi— ^i ^ v ^i...i/i > 

/^)Mi— Mi /n ai "' Q i _ o^ 1 — ^ ('Q9 N 1 

Taking into account the definition of projection operators (1921) and the properties of the X- 
operators (J9~lj) . we obtain 

^ • • ■ ^MiC™ = (93) 

This equation is the basic property of the projection operator: it projects any operator with L 
indices onto the partial wave operator with angular momentum L. 
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For the lowest states, 

= 1, OS = g^, 

_ I ( g ± Q ±- +(] ± Q ± __ Q -L -L ^ f94) 

For higher states, the operator can be calculated using the recurrent expression: 

1 / L 

r\fii...lJ,L — [ „-L r)Mi-A*i-iMi+l-Mi (o^ 

^ ;/!...!/£, r 2 I yfj, i Uj'^v 1 ...Vj- 1 Vj +1 ...i' L W'-'J 



L2 ^=1 



(2L - 1)(2L - 3) 



L 

yA t i/iji'v fe i/ m l - y !yi...!y fc _ 1 iy fe+1 ...iy m _iiy m+ i...i' L 



The product of two X-operators integrated over a solid angle (that is equivalent to the inte- 
gration over internal momenta) depends only on the external momenta and the metric tensor. 
Therefore, it must be proportional to the projection operator. After straightforward calcula- 
tions we obtain 

/§4tl„J* ± )^Lj^) = ^0« V - . (96) 

Let us introduce the positive valued \k\ 2 : 

|g|2_ fc 2 _ [s-(m 1 +m 2 ) 2 ][s-(m 1 -m 2 ) 2 ] ^ 

In the c.m.s. of the reaction, k is the momentum of a particle. In other systems we use this 



definition only in the sense of \k\ = y—k 2 ^; clearly, \k\ 2 is a relativistically invariant positive 
value. If so, equation fl96|) can be written as 



2L 



[—Xf\ (fc^lf) (k^) = ^d!^(-l) L 0^. (98) 

J 4:71 ^ 1 --- flL v\...v L \ ! 2L-\-l Vl—VL K 1 

The tensor part of the numerator of the boson propagator is defined by the projection operator. 
Let us write it as follows: 

F^ = (~l) L 0^ , (99) 
with the definition of the propagator 

K < 100 > 

This definition guarantees that the width of a resonance (calculated using the decay vertices) 
is positive. 
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Figure 2: The acceptance-corrected Dalitz plot for the pp annihilation into 7r 7r°7r° in liquid H2 
and the result of the two solutions. 
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Figure 3: Mass projection of the acceptance-corrected Dalitz plot and angular distributions for 
specific mass slices. The data are taken from3lke pp annihilation into 3tc° in liquid if 2- Solid 
curves correspond to the solution 1 and dashed curves to the solution 2. 
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Figure 4: The acceptance-corrected Dalitz plot for the pp annihilation into 7r°7r 7r° in gaseous 
H 2 and the result of the two solutions. 




Figure 5: Mass projection of the acceptance-corrected Dalitz plot and angular distributions for 
specific mass slices. The data are taken from3fhe pp annihilation into 7r°7r 7r° in gaseous H 2 . 
Solid curves correspond to the solution 1 and dashed curves to the solution 2. 
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Figure 6: The acceptance-corrected Dalitz plot for the pp annihilation into rjir ^ in liquid H2 
and the result of the two solutions. 
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Figure 7: Mass projections of the acceptance-corrected Dalitz plot and angular distributions 
for specific mass slices. The data are taken from the pp annihilation into r)7r°r) in liquid H 2 . 
Solid curves correspond to the solution 1 and dashed curves to the solution 2. 




Figure 9: Mass projections of the acceptance-corrected Dalitz plot and angular distributions 
for specific mass slices. The data are taken from the pp annihilation into r)7r°r] in gaseous H 2 - 
Solid curves correspond to the solution 1 and dashed curves to the solution 2. 
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Figure 10: The acceptance-corrected Dalitz plot for the pp annihilation into ir 7T°ri in liquid H 2 
and the result of the two solutions. 




Figure 11: Mass projections of the acceptance-corrected Dalitz plot and angular distributions 
for specific mass slices. The data are taken from the pp annihilation into 7i 7r°r) in liquid H2. 
Solid curves correspond to the solution 1 and clashed curves to the solution 2. 
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Figure 12: The acceptance-corrected Dalitz plot for the pp annihilation into 7r°7r°?7 in gaseous 
H 2 and two solutions. 
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Figure 13: Mass projections of the acceptance-corrected Dalitz plot and angular distributions 
for specific mass slices. The data are taken from the pp annihilation into 7r°7r°?7 in gaseous H 2 . 
Solid curves correspond to the solution 1 and dashed curves to the solution 2. 
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Figure 14: Solution 1. The contributions of S-wave (two left columns) and D-wave (three 
right columns) to loo moment integrated over t intervals. First line: —0.1 < t< —0.01 GeV 2 , 
second line: -0.2 < t < -0.1 GeV 2 , third line: -0.4 < t < -0.2 GeV 2 and the bottom line: 
-0.4<t<-1.5 GeV 2 . 
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Figure 15: Solution 2. The contributions of S-wave (two left columns) and D-wave (three 
right columns) to Foo moment integrated over t intervals. First line: —0.1 <t< —0.01 GeV 2 , 
second line: -0.2 < t < -0.1 GeV 2 , third line: -0.4 < t < -0.2 GeV 2 and the bottom line: 
-0.4<t<-1.5 GeV 2 . 
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Figure 16: The description of the moments extracted at —0.1 <t< —0.01 GeV 2 (the left two 
columns) and — 0.2 <t<— 0.1 GeV 2 (the right two columns). Dashed curves correspond to the 
solution 1 and full curves to the solution 2. 




Figure 17: The description of the moments extracted at —0.4 <t< —0.2 GeV 2 (two left columns) 
and — 1.5 < t < —0.4 GeV 2 (two right columns). Dashed curves correspond to the solution 1 
and full lines to the solution 2. 
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Figure 18: Pole structure of the 2 ++ -amplitude in the region of the u;a;-threshold: the resonance 
/ 2 (1560): a) solution l(i), b) solution l(ii) 
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Figure 19: From left to right: a) The tttt — > tttt S-wave amplitude squared, b) the amplitude 
phase and c) the Argand diagram for the S-wave amplitude tttt — > tttt. 
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Figure 20: From left to right: The tttt — ► txtx D-wave amplitude squared, the amplitude phase 
and the Argand diagram for the amplitude. 
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Figure 21: Angular distributions for specific mass slices, a) pp — * 3ir° (gas) b) pp — > ■k q t]T] 
(liquid), c,d,e) pp — > 7i rjrj (gas). Solid curves correspond to the solution 2 and dashed curves 
to the solution 2(-) with excluded /o(1300). 
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Figure 22: The description of the moments extracted at —0.1 < t < —0.2 GeV 2 (two left columns) 
and — 0.4 <t<— 0.2 GeV 2 (two right columns). Solid curves correspond to the solution 2 and 
dashed line to the solution 2(-) with excluded /o(1300). 
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Figure 23: Solution 2(-) with / (1300) excluded from the fit. The contributions of S-wave 
to Y 00 moment integrated over different t intervals. First line: t < —0.1 GeV 2 , second line: 
-0.1 < t < -0.2 GeV 2 , third line: -0.2 < t < -0.4 GeV 2 and the bottom line: -1.5 < t < -0.4 
GeV 2 . 
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Figure 24: From left to right: The wir — > txtx G-wave amplitude squared, the amplitude phase 
and the Argand diagram for the amplitude. 
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